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SUMMARY 


The  momentumless  wake  of  a  slender  propeller-driven  model  in  a 
wind  tunnel  at  a  Reynolds  number  of  6.3  x  10^  (based  upon  diameter)  is 
examined  using  the  hot-wire  anemometer.  Measurements  of  the  mean  flow 
field  including  swirl  are  presented,  as  well  as  measurements  of  the  turbu¬ 
lent  intensity  and  Reynolds  stress  witnin  the  first  20  diameters  behind 
the  model.  The  data  are  compared  with  a  previous  experiment  at  a  similar 
Reynolds  number  but  which  used  a  blunt  body  and  a  high-speed  jet  as  a 
propulsion  source.  Significant  differences  are  noted  in  the  development 
of  the  two  wakes  which  are  caused  by  body  bluntness  and  propeller-induced 
swirl.  For  example,  the  turbulent  intensity  is  lower  by  a  factor  of  two 
for  the  present  data  and  the  wake  width  is  somewhat  narrower.  The  effects 
of  swirl  on  the  mean  flow  field  are  important  for  x/D  <  6,  whereas  farther 
downstream  the  swirl  component  primarily  affects  the  turbulent  energy 
balance.  A  theoretical  model  of  the  flow  is  also  presented  which  describes 
the  flow  quite  well,  given  the  measured  conditions  at  one  location.  The 
constants  of  the  turbulence  model  are  determined  experimentally  and  com¬ 
pared  wit*  determinations  made  'or  other  flow  fields.  The  present  values 
of  these  constants  are  comparable,  although  the  strict  universality  of 
these  constants  is  not  conclusive. 
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NOMENCLATURE 


Cross-sectional  area  of  model,  uRg2 
Turbulent  wake  radius 
Normalize  wake  radius.  b*/Rg 
Constants  defined  in  Equations  1  -  4 
Drag  coefficient, drag/q^Ag 
Pressure  coefficient,  (p-pj/q^ 

Torque  coefficient,  torque/q^AgRg 
Maximum  diameter  of  model 
Propeller  diameter 

Turbulent  energy,  j  ( u ' -  +  v'2  +  w'2) 

Curvature  of  normalized  velocity  profile  on  q_ 

Propeller  speed  parameter,  nDp/U^ 

Constant  defined  in  Equation  4 

3/2 

Dissipation  constant,  eA/e 
Eddy  viscosity  constant,  Vy//e  •  A 
Propeller  speed 
Static  pressure 
Pitot  pressure;  p  +  q 
Dynamic  pressure,  j  pu2 
Distance  from  wake  centerline 
Maximum  radius  of  model 
Reynolds  number,  UJJ/v 
Turbulent  Reynolds  number  /iX/v 
Autocorrelation  function 


C<> 


TV, 


u 1  w' 


X,  y,  z 


c 


r, 


A 

A 


v 


VT 

£ 

P 


T 


Subscripts 

o 


ref. 
x,  y 
s 


Veloc1*  cy 

Axial  velocity  difference  from  U^,  maximum 
value  at  given  station 

Swirl  velocity,  maximum  value  at  given  station 
Turbulent  velocity  in  x,  y,  z  directions 
Normalized  RMS  turbulent  velocities  Vo  *  2/U  ,  etc. 

»  CD 

Turbulent  Reynolds  stress 

Distances  measured  frcm  propeller  plane,  x  positive 
downstream,  y  positive  in  North  direction,  z  positive 
in  East  direction 

T  ioulent  dissipation  rate 

Propeller  efficiency 

Longitudinal  integral  scale  of  turbulence 
Longitudinal  microscale  of  turbulence 
Kinematic  viscosity 
Eddy  viscosity 

Normalized  axial  dimension  x/D 
Fluid  density 
Time  lag  for  R ( t ) 

Characteristic  time  of  R(t)  defined  in  Equation  12 

Value  on  centerline  of  wake 

Value  in  freestream 

Reference  or  initial  condition 

Quantity  measured  in  x  or  y  direction 

Swi/1 
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1.  INTRODUCTION 


The  turbulent  wake  created  by  a  self-propelled  body,  such  as  an 
airplane  or  a  submarine,  has  recently  been  a  subject  of  intense  research 
because  of  its  relationship  to  the  question  of  submarine  detection. 

For  example,  in  an  ocean  thermocline  a  submarine  wake  can  sufficiently 
stir  the  fluid  so  as  to  cause  a  sudden  vertical  collapse  of  the  wake  which, 
in  turn,  generates  internal  gravity  waves  (Reference  1).  Unfortunately, 
very  little  is  Known  about  the  so-called  momentumless  turbulent  wake 
either  experimentally  or  theoretically,  so  that  calculations  of  wake-induced 
phenomena  are  unreliable. 

It  is  surprising  that  so  little  experimental  data  exist  even  for 
the  case  of  a  momentumless  wake  in  a  uniform  fluid.  Only  two  studies 
treating  the  self-propelled  wake  have  been  published  -  that  by  Naudascher 
(Reference  2)  and  by  Ginevsxii,  et.  al  (Reference  3).  In  each  of  these 
experiments  the  "body"  was  a  circular  disk  norma1  to  the  oncoming  fluid 
and  the  method  of  thrust  generation  was  a  high-speed  axial  jet.  This 
configuration  does  not  necessarily  provide  a  good  simulation  of  a  full 
scale  flow  field  principally  because  the  model  is  not  slender  or  streamlined 
and  also  because  a  high-speed  jet  is  not  a  common  propulsion  mechanism. 

It  is  well  known  that  a  blunt  body  can  cause  large  scale  fluctuating 
disturbances  to  the  passing  fluid.  An  example  of  such  a  disturbance  is 
the  often  studied  Karman  vortex  street.  Even  at  high  Reynolds  numbers  where 
the  entire  flow  surrounding  the  body  is  turbulent  this  vortex  street 
phenomenon  can  still  be  observed  (Reference  4).  The  turbulent  wake,  on 
the  other  hand,  is  composed  of  smaller  scale  eddies  on  the  order  of  a 
fraction  of  the  wake  radius.  Thus,  for  a  blunt  body,  the  large  scale 
body-induced  disturbances  could  mask  the  development  of  the  turbulent 
wake.  A  streamlined  body,  in  contrast,  creates  much  less  of  a  gross  flow 
disturbance  and  the  wake  develops  in  a  more  realistic  manner. 
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The  experiments  of  Naudascher  and  Ginevskii,  moreover,  add  mass  to 
the  flow  by  virtue  of  the  jet  employed  in  their  test  configuration.  To 
t.he  external  flow  field  this  mass  addition  is  equivalent  to  the  flow  past 
a  point  source  which  has  a  thrust  equal  to  pUjn  (where  m  is  the  source 
strength).  The  effect  of  this  flow  perturbation  on  either  the  force 
balance  or  the  wake  development  is  not  known. 

A  high  speed  jet  supplies  more  energy  to  the  flow  than  is  necessary 
even  though  the  resultant  thrust  is  equal  to  the  body  drag  (i.e.,  it  is 
less  efficient  than  most  other  propulsion  devices).  An  estimate  of  the 
propulsive  efficiency,  o,  can  be  estimated  from  the  relation 
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where  aV  is  the  change  in  fluid  velocity  imparted  by  the  propulsion  device. 
The  excess  energy  is  eventually  converted  to  turbulence  in  regions  of  high 
shear.  Thus,  it  miqht  be  expected  that  the  turbulent  intensities  will 
be  lower  for,  say,  a  propeller-driven  body  than  for  a  jet-driven  body. 

The  use  of  a  propeller  as  a  propulsive  mechanism  is  commonplace  for 
ocean  vessels  because  of  its  relatively  high  efficiency.  An  additional 
complication,  is  that  a  propeller  does  induce  a  swirl  component  within  the 
wake  which  can  have  a  pronounced  effect  on  its  development  such  as  its 
growth  and  the  decay  of  the  turbulent  intensities. 

Because  of  the  deficiencies  noted  in  the  previous  experiments  and 
because  of  the  paucity  of  experimental  date ,  it  was  felt  that  a  more 
realistically  created  momentumless  wake  should  be  experimentally  studied. 
This  report  describes  the  results  of  such  an  experimental  program. 

The  results  of  the  present  experiment  on  the  wake  of  a  slender 
propeller-driven  body  are  presented  in  Section  3.  These  results  are 
compared  with  previous  data  in  an  effort  to  assess  the  differences  between 
the  two  sets  of  data  and  to  determine  if  these  differences  are  caused  by 
body  shape  or  propeller-induced  swirl. 
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A  set  of  turbulent  model  equations  are  also  solved  to  describe  the 
flow  using  the  measured  initial  conditions  at  one  station  in  the  wake.  This 
theory  is  given  in  Section  3.1  and  the  result  compared  with  the  data  in 
subsequent  sections.  A  comparison  of  the  universality  of  the  three 
constants  obtained  from  this  experiment  is  made  with  other  free  turbulent 
flows  in  Section  4. 


2.  TEST  FACILITY  AND  APPARATUS 


2.1  WIND  TUNNEL  AND  MODEL 

This  experiment  was  performed  in  the  TRW  low-speed  wind  tunnel  at 
a  freestream  velocity  of  64  ft/sec.  This  facility  has  a  test  section 
nominally  11  inches  square  by  88  inches  in  length.  The  cross-sectional 
area  of  the  tunnel  does  increase  in  the  flow  direction  to  compensate  for 
wall  boundary  layer  growth  and  result  in  constant  freestream  velocity. 

The  resultant  freestream  velocity  was  uniform  to  within  1%  and  the 
turbulence  level  was  measured  to  be  roughly  0.25%. 

The  model,  shown  in  Figure  1,  is  a  2.00  inch  maximum  diameter 
Rankine  ovoid  with  a  12.5  degree  conical  boat-tail  and  with  an  overall 
length  of  22  inches.  A  three-bladed  propeller,  also  shown  in  Figure  1, 
was  used  to  provide  the  thrust  and  is  constructed  from  brass  with  a 
diameter  of  1.05  inches.  The  pitch-to-diameter  ratio  of  the  propeller 
blades  varies  from  about  0.95  at  the  hub  to  1.5  at  the  blade  tip  and  the 
blades  themselves  are  0.03  inches  thick  and  have  a  total  area  of  33% 
of  the  propeller  disk  area.  The  propeller  was  driven  by  a  high-speed 
air  turbine  capable  of  speeds  up  to  100,000  RPM. 

The  model  was  suspended  by  means  of  a  ball  joint  at  the  end  of  a 
0.75  inch  O.D.  pipe  which  extends  through  the  convergent  section  of  the 
wind  tunnel.  This  pipe  also  serves  as  the  supply  and  exhaust  lines  for  the 
turbine  air  supply.  Four  small  wires  (0.014  inch  diameter)  were  used 
to  position  and  stiffen  the  pipe  aaainst  vibrations  and  were  located  at 
the  test  section  entrance.  A  photograph  of  the  wind  tunnel  installation 
is  shown  in  Figure  2.  Here,  the  flow  is  vertically  downward  with  the 
compass  orientation  noted  in  the  fiqure.  For  future  reference  the  origin 
of  the  coordinate  system  is  located  at  the  plane  of  the  propeller  with  x 
measured  downstream,  y  in  the  northerly  direction  and  z  in  the  easterly 
direction. 

The  Reynolds  number  based  on  the  freestream  conditions  and  the  model 
diameter  is  6.3xl04  for  the  present  experiment.  This  is  some  25%  qreater 
than  the  Revnolds  number  of  the  Naudascher  experiment. 
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2.2  INSTRUMENTATION 

Turbulence  measurements  and  mean  flow  data  were  obtained  using  a 
ThermoSystems  constant  temperature  anemometer  (Model  1051).  The  sensors 
were  single  and  crossed-wire  probes  made  of  0.00015  inch  diameter  tungsten 
wire  with  a  length  of  0.050  inches.  All  wires  had  a  cold  resistance  of 
5-7  ohms  and  were  operated  at  resistance  50%  greater  than  the  cold  value. 

The  output  of  the  anemometer  was  not  linearized  because  of  the 
low  fluctuation  levels  M0%).  Analytic  corrections  were  made  for  the 
mean  flow  measurements  using  an  in-place  wire  calibration  and  King's 
law  -  both  of  which  were  verified  by  static  and  dynamic  calibration  in  an 
open  jet. 

Auxiliary  equipment  used  with  the  anemometer  included  a  Dymec 
Integrating  digital  voltmeter  (1  second  averaging  time)  (Model  2401 B),  a 
Ballantine  True  RMS  meter  (Model  320A),  a  ThermoSystems  Correlator 
(Model  1015C),  a  ThermoSystems  RMS  meter  (Model  1060),  and  a  Saicor 
Correlation  and  Probability  Analyser  (Model  SAI-42A). 

In  addition  to  the  anemometer  data,  mean  flow  information  was 
also  obtained  using  a  Pitot  tube  in  conjunction  with  Pace  differential 
pressure  transducer  (Model  P90D)  and  a  Whittaker  Carrier  Demodulator 
(Model  CD  10).  The  demodulator  output  was  monitored  on  the  abovementioned 
digital  voltmeter.  The  propeller  speed  was  measured  using  a  General  Radio 
Corp.  Strobotac  (Type  1531 -A). 
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3.  RESULTS 


3.1  THEORETICAL  MODEL 

One  important  result  of  the  present  investigatin'  is  the  concurrence 
of  the  present  data  with  one  currently  popular  tur'  ,’1ence  model.  This 
model,  which  is  a  simplified  version  of  that  prr  josed  by  Ng  and  Spaulding 
(Reference  5),  assumes  that  the  shear  stress  can  be  computed  from  an  "eddy 
viscosity"*  such  that  ^  /e  •  n  (originally  attributed  to  Prandtl  and 
Kolmogorov)  and  that  the  rate  of  dissipation  of  turbulent  energy  is 
e  'v  e  '  /A  where  e  is  the  local  turbulent  energy  (=  ^  u^  1 2)  and  a  is  the 
longitudinal  integral  scale  of  the  turbulence.**  The  simplification 
introduced  here  is  that  the  integral  scale  A  is  a  fixed  fraction  of  the 
width  of  the  turbulent  region.  Given  these  relations  along  with  the 
experimental  observation  that  vy  =  Vy(x)  only,  relatively  simple  differential 
equations  for  the  growth  rate  of  the  wake,  the  decay  of  the  mean  field  and 
the  decay  of  the  turbulent  energy  can  be  derived. 

Without  going  into  detail  these  are: 

(Integrated  r2  •  Momentum  Equation) 

buUD  =  constant  (1 ) 


(<£  Momentum  Equation) 


=  -4Cieo 


(2) 


*  The  author  is  aware  of  the  limitations  of  such  a  concept  for  certain 
flows  such  as  the  wall  jet. 

**This  is  strictly  valid  only  for  large  turbulence  Reynolds  number 
r  =  ZlA  .  For  the  present  tests  R,  >  300  which  is  sufficient  to  justify 
the  use  of  the  stated  dissipation  law. 


6 


(Conservation  of  Angular  Momentum) 
b3Vs  =  constant 


(3) 


(integrated  Turbulent  Energy  Equation) 

^  (eot2)  =  2C2e03/2b  +  C 3 (4) 


""Cre  Cj  i  Kjf0-|  (A0/b) 

C2  =  K  •  K£/A0/b 
K  =  vj/e  .  A 

V  t 

3/  2 

K£  =  e-Ao/e^ 

K  =  constant  from  averaging 

e  and  e  over  the  entire  wake  («0.77),  C3  =  constant  from  relating  the 
angulir  momentum  integral  to  the  maximum  swirl  velocity  Vs  and  |f0"| 
is  the  curvature  of  the  axial  velocity  profile  on  the  wake  axis.  In 
the  derivation  of  these  equations  the  effect  of  swirl  on  the  mean  flow 
field  is  neglected  as  is  the  production  of  turbulence  from  the  mean  axial 
flow  field  —  both  of  which  decay  like  f1/7-  relative  to  the  retained  terms. 

As  a  further  simplification,  the  turbulence  production  term 
C3d(Vs2b2)/d£  is  taken  proportional  to  e03/2  b  so  that  the  turbulent 
energy  equation  may  be  written 


^  (e0b2)  «-2C2*  e03/2  b  (5) 

Equations  (1),  (2),  (3)  and  (5)  can  then  be  integrated  in  closed 
form  to  give: 

b/bREF  ■  [(!♦  <2C1+C2*)  ^  -  5REF)] 2Cl/C2*  *' 
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u0/uR[P  *  [(It  (2C  tc2*)  -M  u 

°REF 


4C!/C2* 

n  -2C]7C7  +1 
(*REF'J 


Cq/Cq 


REF 


WsREF 


ST 

C0+  (2CJ+C2*)  - r-— -  U 
°REF 

[d+  (2C!+C2*)  (5 

°REF 


2C1/C2*  +2 

,  n  -2Ci/C2*  +1 
f*REF;j 

3C!/C2* 

t  n  -2C1/C2*  +1' 
**REF'J 


(6) 


Using  the  measured  values  of  =  0.4,  K£  =  0.53,  A0/b  =  0.21,* 
and  I f 0"  |  ~  8.3  so  that  Cj  =  0.767,  C2  =  1.94  and  C2*  =  1.61  along  with 
the  reference  conditions  taken  at  x/D  =6.0  from  the  present  experiment 
gives 

b(U/b(6)  =  i,*  0.213  (5-6)]  °*244 

MO/M6)  =  [1+0.213  (5-6)]  _0-98 

(7) 

&o(0/&o(6)  «  [1+0.213  (5-6)]"0*755 

Vs(0/Vs( 6}  =  [1+0.213  (c-6)j‘°*732 

Plots  of  these  relations  are  given  in  Figures  5,  8,  10  and  13 
along  with  the  experimental  data  and  the  comparison  between  the  two  is 
seen  to  be  excellent  over  the  region  surveyed. 

3.2  EXPERIMENTAL  RESULTS 
3.2.1  Non -Ax 1  symmetry  of  Wake 

Despite  considerable  effort  a  perfectly  round  wake  was  never 
achieved  during  the  course  of  this  experiment.  An  example  if  the  resul¬ 
tant  momentumless  wake  shape  is  given  in  Figure  3  where  a  map  of  the 
turbulent  energy  contours  is  shown  at  x/D  =15.  As  can  be  seen,  the  wake 


*The  values  of  these  constants  computed  from  other  experimentally 
realized  flows  is  given  in  Section  4. 
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is  elongated  in  the  NW-SE  direction.  This  same  shape  was  found  at  all 
survey  stations  despite  changes  in  either  model  position  or  the  model 
angle-of-attack  (by  swiveling  about  the  ball  joint  connection  to  the 
pipe). 


With  the  model  removed  from  the  test  section  no  appreciable  non¬ 
uniformity  could  be  detected  in  the  frees tream  flow  which  could  account 
for  this  behavior.  Furthermore,  attempts  to  correct  the  problem  such  as 
strut-mounting  the  model,  forcing  the  flow  through  a  screen  at  the  test 
section  entrance,  or  lowering  the  model  farther  into  the  test  section 
were  all  unsuccessful . 

The  cause  for  the  asymmetry  is  believed  by  the  author  to  be  a 
small  cross-flow  at  the  bow  of  the  model  which  is  related  to  the  non- 
symmetric  convergent  section  upstream  of  the  test  section. 

Mean  flow  surveys  of  the  wake  (both  hot-wire  and  Pitot)  show  that 
the  model  drag  (fluid  velocity  less  than  freestream)  is  contained  at  the 
two  ends  of  the  elongated  wake  whereas  the  propeller  thrust  (fluid 
velocity  greater  than  freestream)  is  distributed  over  the  midc^e  of  the 
wake.  Wake  surveys  were  generally  taken  along  the  two  lines  i  icated 
in  Figure  3  once  the  center  had  been  found  by  a  preliminary  inarch. 

Entire  maps  of  the  flow  were  made  at  several  stations  to  determine  the 
momentum  defect  and  to  obtain  the  turbulent  map  such  as  shown  in  Figure  3. 

It  should  be  mentioned  that  the  wake  flow  was  extremely  sensitive 
to  the  model  angle-of-attack.  A  pitch  or  yaw  angle  of  0.2  degrees  could 
easily  be  observed  at  x/D  =  8  because  the  propeller  jet  (or  wash)  would 
be  off  to  one  side  of  the  body  wake.  The  final  model  alignment  was 
achieved  by  repeating  surveys  at  x/D  =  8  until  the  most  symmetric  wake 
was  found. 

3.2.2  Determination  of  Propeller  Speed  for  Momentumless  Wake  Condition 

Pitot  pressure  maps  were  obtained  at  x/D  =  6  and  x/D  =  10  fjr 
several  propeller  speeds  in  order  to  determine  when  the  resultant  propel¬ 
ler  thrust  equaled  the  body  drag.  Assuming  that 
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(8) 


Vu»  =  7  W)/q 


these  maps  were  used  to  compute  the  net  drag  in  the  wake  by  manually 
integrating  with  a  planimeter.  The  results  of  these  calculations  are 
shown  in  Figure  4  where  the  wake  drag  is  normalized  by  the  body  drag 
(Cq  ~  0.29  determined  by  a  similar  procedure)  and  the  propeller  speed  is 
given  in  terms  of  the  advance  ratio  J  =  nDp/U^.  All  subsequent  measure¬ 
ments  were  made  at  a  propeller  speed  corresponding  to  J  =  1.58  ±  0.005. 

The  actual  propeller  speed  at  this  condition  is  around  67,000  RPM  which 
was  held  constant  to  within  ±300  RPM  by  regulating  the  turbine  air  supply 
and  monitoring  with  the  Strobotac.  A  tribute  is  due  the  manufacturer  of 
the  air  turbine*  in  that  so  high  a  speed  could  be  easily  maintained  over 
the  course  of  the  measurements  (  150  to  200  hours  or  more).  Further,  one 
can  imagine  the  difficulties  encountered  in  balancing  the  propeller  and 
drive  shaft  at  such  high  speeds  which  was  accomplished  by  a  trial  and 
error  approach. 

3.2.3  The  Mean  Flow  Field 

The  measurements  of  the  wake  radius  obtained  from  the  hot-wire 
surveys  are  shown  in  Figure  5.  The  wake  radius,  b*,  is  defined  here  as 
the  location  where  the  axial  turbulent  intensity  dropped  to  one-half  of 
the  centerline  value.  Data  points  are  shown  for  two  different  probes  and 
the  upper  set  of  points  represent  the  x-axis  measurement  and  the  lower  set 
the  y-axis  measurement  (see  Figure  3  for  orientation).  An  equivalent  or 
mean  wake  radius  defined  as 

b*/RB  =  [bx*.by*/Rg2]i/2 

is  indicated  by  the  solid  circles.  Shown  also  is  the  wake  radius  measured 
by  Naudascher  downstream  of  a  jet-driven  circular  disk.  The  present 
results  are  somewhat  smaller  than  the  disk  wake  results  and  appear  to  reach 
similarity  (b*«  x«)  sooner.  The  growth  of  the  present  wake  is  seen  to  be 


*D0TC0,  Inc.,  Hicksville,  Ohio 
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in  excellent  agreement  with  the  theoretical  curve  indicated  by  the  heavy 
dashed  line.  The  growth  rate  over  the  region  6  <  x/D  <  2C  is  approximately 
x  .  Much  farther  downstream  the  theory  indicates  the  growth  rate  will 
decrease  slightly  to  x°‘"44. 

The  wake  velocity  profiles  Ug/U^  for  x/D  >_  6  are  shown  in  Figure  6 
for  the  x-axis  surveys.  Velocity  profiles  for  x/D  =  2,  3  and  4  are  shown  in 
Figure  7  where  the  near-field  effects  of  the  propel  'r  dominate  (such  as 
swirl-induced  pressure  and  propeller  thrust  distribu  on).  Note  that  close 
to  the  body  that  the  propeller  acts  like  an  annular  jet  and  that  this  near 
field  effect  dies  out  downstream  of  x/D  ~  6. 

The  maximum  velocity  excess  measurements  Up/U  are  shown  in  Figure  8 
where  measurements  other  tnan  those  taken  from  Figure  6  are  presented. 

The  three  determinations  of  Up/U^  are  quite  consistent  with  one  another. 
Indicated  on  this  figure  is  the  near-field  effect  which  causes  the  centerline 
velocity  to  be  less  than  the  peak  velocity  for  x/D  ■'  b.  The  centerline 
velocity  distribution  of  Naudascher  is  also  shown  which  is  initially  much 
higher  than  the  present  experiment  [UD  (x  *  0)/Uqo«  3.64  compared  to  a  value 
of  1.3  or  so  here)  and  exhibits  a  distinctly  different  decay  rate  over  the 
first  20  body  diameters.  The  theoretical  prediction  of  Equation  (7)  is 
seen  to  be  in  excellent  agreement  with  the  data. 

Measurements  of  the  swirl  velocity  (obtained  from  the  crossed-wire 
anemometer  along  the  z-axis)  are  shown  in  Figure  9,  Here,  an  attempt  is 
made  to  display  the  similarity  of  the  profiles  by  normalizing  the  swirl 
velocity  by  the  local  maximum  value  and  normalizing  the  radial  distance  by 
the  turbulent  wake  radius.  The  inner  20  percent  of  the  wake  is  seen  to  be 
in  solid  body  rotation  (v$  ^  r)  and  the  maximum  swirl  velocity  occurs  at 
r/b  ~  0.25  to  0.3.  The  solid  line  faired  through  the  data  has  been  used  to 
compute  the  radial  swirl-induced  pressure  distribution,  the  angular 
momentum  and  the  flux  of  swirl  energy. 


For  the  reduction  of  the  Pitot  tube  data  to  velocity  the  local  static 
pressure  was  computed  using  the  swirl  velocity  data  and  the  turbulent 
intensity  data.  Direct  measurement  of  the  static  pressure  was  not  made 
as  it  was  quite  small  (^0.01  q  ). 

*oo 
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The  maximum  swirl  velocity  measurements  are  shown  in  Figure  10 
along  with  the  theoretical  line  corresponding  to  V$b3  =  constant  (i.e., 
conservation  of  angular  tncmentum).  It  should  be  noted  that  the  swirl 
velocities  are  comparable  in  magnitude  to  the  axial  flow  velocities. 

The  swirl -induced  pressure  distribution  calculated  from 

oo 

P.  -  P<r>  *  »./  %  00) 

r 

and  using  the  faired  line  from  Figure  9  is  shown  in  Figure  11.  The  radial 
swirl-induced  pressure  field  is  normalized  by  the  axis  pressure  which  is 
calculated  to  be  CPq  SW)r^  =  “149  (Vs/Uj  from  Equation  (10).  The 
swirl  induced  pressure  deficit  from  freestream  is  seen  from  Figure  11 

to  be  signiricant  only  in  the  core  of  the  wake  (r/b  <  0.5)^ The  turbulence- 

induced  pressure  deficit  in  the  wake  given  by  CpTURB  a  -2  u'2  is  slightly 
less  on  centerline  than  the  swirl-induced  pressure  deficit. 

The  torque  exerted  by  the  propeller  on  the  fluid  has  also  been 
calculated  and  the  torque  coefficient  defined  as 


r  _  Torque 


(11) 


has  a  mean  value  of  0.075  ±  0.015.  Incidentally,  the  overall  efficiency 
of  the  propeller  defined  as  the  thrust  power  'drag  .  U  )  divided  by  total 
input  power  to  the  propeller  (Torque  .  shaft  speed)  or  npROp  =  ^  (CD/CT-J) 
is  0.41.  This  is  somewhat  lower  than  expected  but  not  unreasonable  in 
light  of  the  large  Reynolds  number  difference  [0(103)]  between  the  model 
and  full  scale  propellers. 

Lastly,  the  rotation  of  the  wake  "core"  has  been  computed  by  integrating 
the  measured  angular  swirl  velocity  along  the  wake.  As  the  flow  proceeds 
from  3  to  20  diameters  behind  the  body,  the  inner  part  of  the  wake  revolves 
nearly  0.8  revolutions.  Closer  to  the  propeller,  of  course,  the  swirl  is 
more  intense  causing  perhaps  2  to  6  revolutions  by  x/D  =  3.  This  swirling 
motion  could  be  very  important  in  mixing  up  a  wake  inasmuch  as  this  is  a 
gross  convective- type  motion  as  opposed  to  the  smaller  scale  turbulent 
mixing. 
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3.2.4  The  Turbulent  Field  of  the  Wake 

Tne  axial  turbulent  intensity  (Vu^/U  =  u)  distributions  are  shown 

oo 

in  Figure  12  where  the  centerline  intensity  and  the  wake  radius  have  been 
used  to  normalize  the  intensity  and  radial  position  respectively.  The 
maanitude  of  the  centerline  intensity  uQ  is  shown  in  Figure  13  where, 
again,  other  measurements  are  also  included.  Note  that  the  present  results 
are  a  factor  of  2  or  so  lower  than  observed  by  Naudascher  for  be  jet-driven 
disk  wake.  This  is  a  consequence  of  the  fact  that  a  large  propeller  is  a 
more  efficient  propulsion  device  than  a  small  high-speed  jet.  The  decay 
rate  of  the  turbulent  intensity  is  nearly  the  same  for  both  experiments, 
however.  The  theoretical  prediction  given  by  Equation  (7)  is  seen  to  predict 
the  decay  rate  quite  well  although  it  does  underestimate  the  data  by  5% 
or  so. 


The  transverse  components  of  turbulence  v  and  w  were  also  measured 
by  means  of  the  crossed  wire  and  a  typical  profile  is  shown  in  Figure  14 
taken  at  x/D  *  8  along  with  a  line  representing  the  axial  fluctuations,  u. 
Although  there  are  noticeable  differences  between  u  and  v  or  w,  the  overall 
impression  is  that  the  three  intensities  are  equal  over  .he  bulk  of  the  wake. 
This  fact  will  be  used  to  relate  the  turbulent  energy  e  to  the  axial 
intensity  (i.e.,  e1/2  «  /37?u). 

Measurements  of  the  Reynolds  stress  u'v'  are  shown  in  Figure  15  for 
transverses  taken  along  the  x-axis.  The  asymmetry  of  these  profiles  corresponds 
to  the  observed  asymmetry  of  the  mean  velocity  profiles  (Figure  6). 

In  addition  to  determining  the  turbulent  intensities  and  Reynolds 
stresses  the  anemometer  signal  was  also  connected  to  the  Correlation  and 
Probability  analyzer  for  a  more  thorough  statistical  analysis.  Auto¬ 
correlation  functions  of  the  axial  turbulence  component  were  taken  at  a 
number  of  stations  along  the  wake  axis  and  also  outward  along  a  radius.  A 
typical  example  of  the  u'  autocorrelation  function  R( T )  is  shown  in 


~t/ ' 


Figure  16  for  x/D  =  10.  Note  that  R(x)  is  closely  represented  by  e 

and  that  no  periodic  behavior  is  observed.* 

* 

h!!\ppn  th!  p6ak  w-1  v!loc1>  location  a  weak  periodic  signal  could 

fLonf  /ith?  J|Jt°corj[elat''or-  funct1on  corresponding  to  the  propeller 
frequency  (1.1  KHz).  By  x/D  =  3  this  behavior  could  no  longer  be  detected-- 
it  presumably  being  destroyed  by  the  turbulence. 
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The  coirmon  definition  of  the  integral  scale  of  turbulence  is 
(accepting  Taylor's  hypothesis) 

oo 

Ao  5  r(t'  <*t  (12) 

when  R ( t )  =  e"T//l  this  becomes 

VD  =  V  (]3) 

The  measurements  of  AQ/D  determined  in  this  manner  are  shown  in  Figure  17 
for  the  wake  centerline.  Shown  also  are  the  data  of  Naudascher  for  the 
jet-driven  disk  which  are  somewhat  greater  than  found  here.  The  measure¬ 
ments  of  A  have  also  been  normalized  by  the  local  wake  radius  and  the  fact 
that^these  two  are  proportional  is  shown  bv  the  solid  like  which  indicates 
A0/b  ~  O.ti.  This  is  still  somewhat  smaller  than  found  by  Naudascher 
who  reports  A  /b*  «  0.28  for  x/D  £  20  but  decreasing  to  0.21  by  x/D  %  50. 
The  fact  that  the  present  experiment  shows  A/b  =  constant  has  been  used 
in  the  theoretical  analysis  previously  given. 

Incidentally,  the  turbulent  integral  scale  size  was  also  measured 
for  the  drag  wake  and  a  value  ol  A/b  ~  0.25  ±  0.02  was  determined. 

Acnss  the  momentumless  wake  the  turbulent  scale  size  wai  observed 
to  increase  somewhat  toward  the  edge  of  the  wake.  This  increase  is 
undoubtedly  effected  by  intermi ttency  (which  was  not  measured)  and  is 
consistent  with  the  assumption  that  \>y  is  constant  across  the  wake  (e.g., 
since  vy  *  /e  A~*e^  aq  then  a/aq  =  [e/eQ]'1/2). 

The  probability  density  of  axial  fluctuations  and  transverse 
fluctuations  were  also  measured  and  found  to  be  very  close  to  a  normal 
distnbut1on--at  least  on  the  wake  axis.  This  result  was  found  at  all 
stations  even  as  close  as  x/D  =  2.  Probability  density  measurements  away 
from  the  wake  axis  were  skewed  which  reflects  the  intermittent  nature  of 
the  flow  near  the  wake  edge. 


14 


An  attempt  was  also  made  to  measure  the  microscale  of  turbulence 
by  measuring  the  quantity  ( 3u ' / yt ) ^  and  using  the  formula  given  by  Hinze 
(Reference  6) 

X2=U2^/(ff)2  (14) 

These  measurements  were,  however,  not  repeatable  from  day-to-day  and 
moreover  appeared  somewhat  large.  This  is  probably  caused  by  the  fact 
that  the  length  of  the  wire  was  +he  same  order  as  microscale.  A  wire 
length  correction  procedure  outlined  by  Wygnanski  and  Fiedler  (Reference  7) 
was  not  used  because  of  the  unrepeatability  of  the  \  measurements 
(30  percent  variation). 

3-3  DETERMINATION  OF  EDDY  VISCOSITY 

The  mean  velocity  profiles  and  the  Reynolds  stress  measurements 
(Figures  6  and  15,  respectively;  were  used  to  compute  the  turbulent  or 
"eddy"  viscosity  through 


u-  u . 
i  i 
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This  procedure  is  admittedly  inaccurate  by  virtue  of  the  fact  that 
experimental  data  must  be  differentiated.  Nevertheless,  its  evaluation 
was  felt  important  enough  to  warrant  the  effort  especially  in  light  of  the 
fact  that  Naudascher's  data  indicate  that  the  eddy  viscosity  is  roughly 
constant  throughout  the  wake  and  which  is  at  odds  with  recent  scaling  laws. 

The  radial  variation  of  the  eddy  viscosity  is  shown  for  one  station 
(x/D  =  6)  in  Figure  18  where,  except  near  the  edge,  it  appears  to  be 
relatively  constant.  The  axial  variation  computed  for  each  profile  at  the 
point  of  maximum  slope  and  maximum  uV  for  each  side  of  the  wake  is  shown 
in  Figure  20  normalized  by  the  laminar  viscosity  v.  The  vertical  bars 
indicate  the  inaccuracy  of  the  calculation.  Note  that  the  turbulent  viscosity 
decreases  by  a  factor  of  2  to  2-1/2  over  the  region  3  <  x/D  <  15,  which  is 
in  contrast  to  the  blunt  body  jet-propelled  result.  Shown  also  is  the 
Prandtl -Kolmogorov  eddy  viscosity  relation  vT  =  Ky  v€a.  For  a  value  of 

K  ~  °*4  this  relation  provides  an  excellent  approximation  to  the 
experimental  data. 
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3.4  ENERGY  BALANCE  IN  WAKE 


The  integrated  flux  of  various  forms  of  energy  (turbulent,  swirl 
and  mean  flow)  along  the  wake  have  also  been  computed.  These  results  are 
shown  in  Figure  20  where  the  various  fluxes  have  been  normalized  by  the 
model  drag  power  (drag  •  U^).  For  example,  the  turbulent  energy  flux  is 
given  by 


*turb  =  /  u  e  do  <16) 

wake 

and  ihe  integral  is  evaluated  by  manual  integration  of  the  contour  map 
shown  in  Figure  3.  Similar  techniques  were  used  to  calculate  the  swirl 
energy  flux  E$wir1  =  w^£  u  (vg2/2)  da  and  mean  flow  energy  flux 

^mean  =  wake  u  ^ud  do‘  Note  ener9y  fluxes  are  quite 

small  compared  to  the  model  drag  power--less  than  10  percent  downstream 
of  x/D  =  4.  Further,  the  flux  of  turbulent  energy  is  far  greater  than 
either  the  swirl  energy  flux  or  mean  flow  energy  flux.  The  decay  rates  of 
the  turbulent  energy  flow  and  the  swirl  energy  flux  are  comparable  whereas 
the  mean  flow  energy  decays  somewhat  more  rapidly  which  is  in  agreement 
with  the  theoretical  model. 

Shown  for  comparative  purposes  is  the  amount  of  turbulent  energy 
deposited  into  the  wake  from  the  boundary  layers  which  was  computed  using 
the  flat  plate  data  of  Klebanoff  (Reference  8  and  taken  from  Reference  9) 
and  using  Schlichting  (Reference  10)  to  estimate  the  boundary  layer 
thickness  and  friction  darg.  This  estimate  shows  that  the  bulk  of  the 
turbulent  energy  is  generated  downstream  of  the  propeller  by  the  axial 
and  swirl  velocity  gradients  rather  than  simply  being  convected  through  the 
propeller  from  the  boundary  layers. 

The  dissipation  rate  of  turbulent  energy  can  be  computed  by  assuming 
that  the  loss  of  swirl  and  mean  flow  energies  is  effectively  a  source  or 
production  term  for  the  turbulent  energy.  Assuming  that  the  dissipation  rate 
of  turbulent  energy  on  centerline  is  proportional  to  eQ3/2/A0,  the  dissipation 
constant  is  found  to  be  0.53  ±  0.02  which  is  in  fair  agreement  with  the 
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value  determined  from  other  experimental  flows--see  Section  4  for  a  more 
detailed  explanation.  The  effective  dissipation  constant  used  in  the 
theoretical  model  is  0.44  which  was  calculated  from  the  measured  turbulent 
energy  flux  curve  shown  in  Figure  21. 
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4.  COMPARISON  OF  EVALUATED  TURBULENT  CONSTANTS 
WITH  OTHER  EXPERIMENTAL  RESULTS 


In  the  formulation  of  a  set  of  turbulence  model  equations  such  as 
that  of  Ng  and  Spaulding,  or  others,  one  has  the  hope  that  the  few  arbitrary 
constants  used  to  characterize  the  turbulence  might  be  universal  in  some 
sense.  In  other  words,  these  constants  should  be  independent  of  tne 
particular  flow  field  under  investigation.  In  the  present  experiment  three 
such  constants  were  evaluated  for  the  case  of  the  momentumless  wake  with 
swirl.  These  are  K^,  the  eddy  viscosity  constant;  ,  the  dissipation 
constant;  and  AQ/b* ,  the  tu-bulent  integral  scale  size  normalized  by  a 
characteristic  dimension  of  the  extent  of  turbulent  flow,  and  it  is  of 
interest  to  see  how  these  values  compare  with  other  experimentally 
realized  flows. 


For  this  comparison  five  other  flows  have  been  chosen  —  two  planar 
flows  (References  11,  12),  two  axi symmetric  flows  (References  7,  13,  14), 
and  homogeneous  grid  turbulence  (References  14,  16)*  The  evaluation  of 
the  five  turbulent  quantities  e,  A,  vp  e  and  b*  have  all  been  evaluated 
from  the  data  where  possible  or  estimated  if  necessary.  The  results  of 
this  comparison  are  shown  in  Table  1.  Also  shown  in  this  table  is  the 
turbulent  Reynolds  number  which  should  be  larger  than  30  or  so  to  insure 

fbat  any  laminar  phenomena  are  indeed  small  compared  to  the  turbulent 
phenomena. 


A  cursory  examination  of  Table  1  indicates  that  complete  universality 
of  the  three  turbulent  constants  has  not  been  achieved.  A  factor  of  3  or  4 
variation  in  is  noted  whereas  the  variation  in  and  A/b*  is  within  a 
factor  of  2.  The  constants  derived  from  the  present  experiment  look 
reasonable  compared  to  other  flows  although  is  somewhat  higher  than 
average  and  A/b*  is  somewhat  lower.  Whether  these  small  differences  are 
caused  by  the  particular  flow  field  or  by  experimental  inaccuracies  is  unknown. 


* 

Toue-aU^h0r  t0  acknowledge  the  assistance  of  Dr. 
IKW  in  locating  these  references  and  the  calculation  of  the 
from  the  unpublished  data. 


I .  E .  A1 ber  of 
parameters 
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5.  CONCLUSIONS 


Based  upon  the  present  experimental  results  for  the  wake  of  a 
slender  propeller-driven  body  the  following  conclusions  are  drawn: 

o  The  wake  development  is  significantly  different  from  previously 
self-propelled  experimental  results  obtained  for  a  blunt  body 
driven  by  a  high  speed  jet.  For  example,  in  the  present  experiment 
tne  wake  is  narrower  and  the  turbulent  intensities  are  lower  by  a 
factor  of  two.  Further,  the  integral  scale  size  is  25%  lower 
and  the  eddy  viscosity  is  observed  to  decrease  with  distance  aft. 

•  The  effect  of  propeller  swirl  on  the  mean  axial  flow  field  is 
confined  to  the  region  x/D  *  6.  Downstream  of  this,  the  swirl 
effects  the  flow  field  by  acting  as  a  source  of  turbulent 
energy.  The  swirl  velocities  were  observed  to  be  larger  than 
the  axial  velocity  deviations.  The  "core"  of  the  wake  (r  <  0  25b) 
was  found  to  be  in  solid  body  rotation. 


From  measurements  of  the  eddy  viscosity,  integral  scale  size 
and  turbulent  dissipation  rate  a  simple  turbulence  model  was 
found  which  quite  accurately  described  the  wake  development  usinq 
measured  initial  conditions  at  x/D  =  6.  This  turbulence  model 
would  predict  that  far  downstream  in  the  wake  that  the  wake 
radius  would  increase  as  x°-24\  the  turbulent  intensity  would 
ecrease  as  x  •  and  the  axial  and  swirl  velocities  would 
decrease  as  x-°-98  and  x-0-78,  respectively. 


•  The  turbulence  constants  determined  in  this  experiment  compare 
adequately  with  determinations  made  in  other  flow  fields  jets 
wakes,  and  shear  layers).  While  it  cannot  be  said  these 
constants  are  universal,  the  variation  from  flow  is  tolerable 
(factor  of  2  or  3). 
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6.  RECOMMENDATIONS  FOR  FUTURE  WORK 


Even  with  the  experimental  results  presented  in  this  report,  the 
understanding  of  the  development  of  a  full  scale  momentumless  wake  is  still 
not  complete.  For  example,  the  degree  of  swirl  is  not  simulated  exactly 
here  and  the  effects  of  stratification  or  measurements  of  the  turbulent 
mixing  were  not  treated  at  all.  It  is  recommended  therefore  that  the 
following  experiments  be  pursued: 

•  Examine  the  effect  of  different  diameter  propellers  and 
propellers  with  differently  pitched  blades.  These  two 
effects  directly  influence  the  axial  and  swirl  velocity 
fields  and  could,  perhaps,  substantially  change  their 
relative  contributions  to  the  decay  of  the  turbulence. 

•  Measurement  of  the  swirl  flow  Reynolds  stress  is  also  of 
importance.  This  effectively  gives  the  rate  at  which  the 
swirl  component  is  converted  to  turbulence.  Very  little 
is  known  of  how  to  model  this  particular  Reynolds  stress 
and  its  determination  could  impact  other  fields  such  as 
aircraft  wake  vortices  or  swirling  jets. 

•  The  initial  degree-of-mixi ng  in  one  wake  could  be 
determined  by  creating  a  gradient  of  a  passive  scalar  in 
the  tunnel  freestream.  Measurement  of  the  gradient  within 
the  wake  is  a  direct  measure  of  how  well  mixed  the  wake 
flow  is.  This  would  have  direct  impact  on  calculations 
which  account  for  the  gravitational  restoring  forces 

for  flows  within  a  stratified  fluid. 
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MOMENTUMLESS  WAKE  EXPERIMENT 
WIND  TUNNEL 


SIZE.  1'  IN.  x 
SPEED  64  FT ''SEC 
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Figure  2.  Photograph  of  Wind  Tunnel  Installation 


Example  of  Observed  Wake  Asymmetry 
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Figure  5 
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AXIAL  POSITION,  x/D 


The  Turbulent  Wake  Radius 
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irl  VElocity  Radial  Prof 


Figure  11.  The  Swirl -Induced  Radial  Pressure  Distribution 
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AXIAL  POSITION,  x'D 


Figure  13.  The  Decay  of  the  Turbulent  Intens 


Figure  14.  Transverse  Turbulence  Profiles 
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AXIAL  POSITION,  x/D 


Figure  17.  Integral  Scale  of  rurbulence  on  the  Wake  Axis 
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INTEGRATED  ENERGY  FLUX/MODEL  DRAG  POWER 


Figure  20.  Energy  Flux  Measurements  in  the  Momentumless  Wake 
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